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Abstract 



It is well-known that the dominant late time behavior of waves propagating 

in . 

on a Schwarzschild spacetime is a power-law tail; tails for other spacetimes 

O 

, have also been studied. This paper presents a systematic treatment of the tail 

in . 

Q\ , phenomenon for a broad class of models via a Green's function formalism and 

q-i establishes the following, (i) The tail is governed by a cut of the frequency 

(3JT), Green's function G(co) along the — Im co axis, generalizing the Schwarzschild 



result, (ii) The uj dependence of the cut is determined by the asymptotic 
but not the local structure of space. In particular it is independent of the 
presence of a horizon, and has the same form for the case of a star as well, (iii) 
Depending on the spatial asymptotics, the late time decay is not necessarily a 
power law in time. The Schwarzschild case with a power-law tail is exceptional 
among the class of the potentials having a logarithmic spatial dependence, 
(iv) Both the amplitude and the time dependence of the tail for a broad class 
of models are obtained analytically, (v) The analytical results are in perfect 
agreement with numerical calculations. 
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I. INTRODUCTION 



In this paper we are concerned with the systematics of the propagation of linearized 
waves on curved background, in particular, the late time behavior ffl. The problem is of 
interest for understanding wave phenomena in relativity [0-0], and possibly in other areas 
of physics. The propagation of linearized gravitational, electromagnetic and scalar waves is 
often modeled by the Klein-Gordon (KG) equation with an effective potential 

D<f>{x, t) = [d 2 t -d 2 x + V{x)] <f>(x, t) = 0, (1.1) 

which includes the Zerilli and Regge- Wheeler equations as special examples [0]. The effective 
potential V(x) describes the scattering of <fi by the background curvature. We are here 
interested in the evolution of in the presence of this scattering. 

On flat space (V = 0), waves travel along the light cone: <f)(x, t) = 4> + (x — t) + <p_ [x + 1), 
but the situation is much more interesting in the presence of a nontrivial background. Fig. 1 
shows the result of a typical numerical experiment evolving (l.l)on0<x<oo (so that 
x may be thought of as a radial variable, see below), with the boundary condition 0(x = 
0, t) = 0, and the condition of outgoing waves for x — > oo. An initial gaussian centered 
at y Q with initial = is allowed to evolve, and is observed at X\ > y Q . The figure shows 
4>(xi, t) versus t. Our interest in this paper is the behavior at large t, but it is useful to first 
place the time dependence in an overall context. 

The generic time dependence of waves observed at a fixed spatial location x\ consists of 
three components, as illustrated in Fig. 1. (i) The prompt contribution (peaks A and C in 
Fig. la) depends strongly on the initial conditions. The first peak A is due to the initial 
pulse propagating from y Q directly to the right, while the later peak C is due to the initial 
pulse propagating from y a to the left, being reflected and inverted at the origin, and then 
propagating to the observation point x\. (ii) Each of these pulses is followed by a trailing 
signal (respectively B and D) caused by the potential V (in the sense that B and D vanish 
when V = 0). The latter pulse D has traversed regions of larger V, and as a result is both 
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more pronounced in amplitude and also more compressed in time. On a longer time scale 
(Fig. lb), it is seen that scattering by the potential (e.g. peak D) merges into a train of 
quasinormal ringing, decaying in an exponential manner, (hi) Finally, there is a tail, which 
decays approximately as a power in t (a straight line in this log-log plot), and it dominates 
at late times. 

The prompt contribution is most intuitive, being the obvious counterpart of the light- 
cone propagation in the V = case, and is in any event expected in the semi-classical limit. 
The ringing is due to a superposition of quasinormal modes [which are solutions to (1.1) with 
the time dependence e~ luJt for some complex u\. These are particularly interesting in that 
the time dependence (both of the quasinormal modes and of the trailing signals B and D) 
observed far away carries information about the structure of the intervening space. In the 
mathematical sense, this is because the quasinormal mode frequencies u are determined by 
V(x). In a more physical sense, it is useful to consider the analogy of electromagnetic waves 
emitted by a source in an optical cavity; the cavity, in providing a nontrivial environment, 
is analogous to the presence of a potential V(x) in (1.1). [In fact, a one-dimensional scalar 
model of electromagnetic waves would be described by a wave equation with a nontrivial 
dielectric constant distribution n 2 (x) ||; such a wave equation would be very similar to 
(1.1) and in fact can be transformed into it [f|.] For a broad band source in a laser cavity 
of length L, the spectrum of electromagnetic waves observed far away would have frequency 
components uo m jnc/L, where j is an integer. In this sense, these quasinormal modes carry 
information about the cavity. In much the same way, the ringing component of the solution 
in (1.1), if observed, would carry significant information about the background curvature of 
the intervening space. For this reason, there have been extensive numerical studies of various 
models. In addition, one particularly intriguing possibility is that the sum of quasinormal 
modes may be complete in a certain sense, as the normal modes of a conservative system 
are complete. If this were the case, one would have a discrete spectral representation of the 
dynamics, with obvious conceptual and computational advantages. The completeness has 



been demonstrated in the analogous case of the wave equation ||10|| , and its generalization 



to the KG equation (1.1) will be reported elsewhere [[□] . 

This paper develops a general treatment for the late time behavior of a field satisfying 
(1.1) with a broad class of time-independent potentials V(x) and the outgoing wave boundary 
condition. Such a system could describe linearized waves evolving in stationary spacetimes 
with or without horizons, and in vacuum or nonvacuum spacetimes. For example, for a 
general static spherically symmetric spacetime, 



ds 2 = g tt (r)dr + g rr (r)dr 2 + r\d6 2 + sin 2 6dip 



;i.2) 



a Klein-Gordon scalar field $ can be expressed as 



1 



$ = ^2-<f>i m (x,t)Yi m (9,(p). 



;i.3) 



Im 



The evolution of 4>i m (x,t) is given by (1.1) with the effective potential 



V(x) = -g tt 
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where x is related to the circumferential radius r by 



x 



9tt 



;i.5) 



and extends over a half-line if the metric is non-singular, but extends over the full line when 
there is a horizon. 

The special case of linearized perturbations of a black hole background has received much 
attention p|-[7|,p^,p!3[1 . In the simplest case of the Schwarzschild spacetime 



ds 2 



1 )dt 2 + 1 dr 2 + r 2 (dd 2 + sin' 

r J V r / 



;i.e) 



the transformation (1.5) leads to 



x = r + 2M log 



/ r 



V2M 



1 



;i.t) 



For a Klein-Gordon scalar field, V is defined by (1.4) and takes the explicit form 

2M" 



V(x) 



1 - 



1.8) 



with S = 0. Electromagnetic and gravitational waves also obey (1.1) and (1.8), but with 
S = 1,2 respectively 0. The boundary condition is outgoing waves at r = 2M (x — > — oo, 
into the black hole) and as r — ► oo (a; — ► oo). With the effective potential (1.8), it is 
well known that the tail decays as £~( 2 '+ 3 ) independent of the spin S of the field For 
Reissner-Nordstrom black holes, the evolution equation can again be put in form of (1.1), 
but with an effective potential slightly different from (1.8). The tail is again found to decay 
as £~( 2i + 3 ) ||. Although a similar tail should exist for Kerr black holes, we are not aware of 
the corresponding analysis in that case. Most of these previous analyses (except Ref. 3) are 
specific to a particular form of the potential. 

Several recent developments make a more general study of the late time tail phenomenon 
of gravitational systems particularly interesting. As studies of the Schwarzschild spacetime 
show that the tail comes from scattering at large radius [|J, it has been suggested || that 
a power-law tail would develop even when there is no horizon in the background, implying 
that such tails should be present in perturbations of stars, or after the collapse of a massless 
field which does not result in black hole formation. In Ref. 6, the late time behavior of 
scalar waves evolving in their own gravitational field or in gravitational fields generated by 
other scalar field sources was studied numerically, and power-law tails have been reported 
(though with exponents different from the Schwarzchild case) . These interesting results call 
for a systematic analysis of the tail phenomenon in (i) nonvacuum, (ii) time-dependent, and 
(iii) nonlinear spacetimes, for which the present work will lay a useful foundation. 

In Section II, we formulate the problem in terms of the Green's function, leading to a 
precise definition of the prompt, quasinormal mode and tail components in terms of con- 
tributions from different parts of the complex frequency plane. In Section III, numerical 
experiments are presented. The analytic treatment is then given in Section IV (potentials 
decreasing faster than x~ 2 ), Section V (inverse square potentials) and Section VI (composite 
potentials consisting of a centrifugal barrier plus a subsidiary term). The last case, though 
technically most complicated, is also the most interesting physically, since it includes, for 
example, the potential in (1.8). Section VII contains a discussion and the conclusion. 
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The highlights of this paper are as follows. Firstly, the analytic results agree perfectly 
with the numerical experiments, including both the time dependence and the magnitude; 
this agreement indicates that the phenomenon of the late time tail has now been quite fully 
understood for a large class of models. Secondly, our work extends, places in context, and 
provides insight into the known results for the Schwarzschild case. For the Schwarzschild 
spacetime, the late time tail has previously been related to either reflection from the asymp- 
totic region of the potential [Q, or to the Green's function having a branch cut along the 
negative imaginary axis in the complex frequency plane In this paper, we demon- 

strate that the cut is a general feature, if V(x) asymptotically (x — ► oo) tends to zero 
less rapidly than an exponential, except for the particular case of the "centrifugal barrier" 
V(x) = 1(1 + l)/x 2 , I = 1,2.... Furthermore, we determine the strength of the cut in terms 
of the asymptotic form of the potential for a broad class of potentials. These results then 
relate the two different explanations of the late time tail in the literature. 

Some specific late time behavior we discovered in this paper are of particular interest. 
Potentials going as a centrifugal barrier of angular momentum I (I is an integer) plus ~ 
^""(logx)^ ((3 = 0,1) are of interest, because these include the Schwarzschild case (a = 
3,(3 = 1). For such potentials, the generic late time behavior is ~ *-(«+«) (tog t)". Thus, for 
(3 = 1, the late time tail is not a simple power law, but exhibits an additional logt factor. 
However, when a is an odd integer less than 21 + 3, the generic leading behavior vanishes. 
For (3 = 0, the next leading term is expected to be £-( 2 '+ 2q - 2 ) ; while for (3 = 1, the next 
leading term is t^( 2 '+ a ) without a logt factor. Interestingly, the well-known Schwarzschild 
spacetime belongs to such exceptions. 

As a collorary, we establish that the form of the late time tail is independent of the 
potential at finite x. In particular, whether the background metric describes a star (so that 
x is itself the radial variable < x < oo), or a black hole [so that x is the variable defined 
in (1.5), — oo < x < oo], the late time tail is the same as long as the asymptotic potential 
V(x — > oo) is the same, thus confirming the suggestion of Ref. 5. 

As a by-product, we point out that numerical calculations are sometimes plagued by a 



"ghost" potential arising out of spatial discretization. The effects of this "ghost" potential 
may in some cases mask the genuine tail behavior — which is often numerically very small 
(see Fig. lb). This problem could have led to errors, or at least ambiguities, in some 
numerical studies, in particular numerical studies of tail phenomena. This subtlety with 
numerical calculations is discussed in Appendix A. 



II. FORMALISM 



The time evolution of a field 4>(x,t) described by (1.1) with given initial data can be 



written as 



<P(x, t)= J dy G(x, y- 1) dt<P(y, 0) + J dy d t G(x, y; t) <j>(y, 0) (2.1) 
for t > 0, where the retarded Green's function G is defined by 

DG(x,y;t) = 5(t)5(x-y) , (2.2) 
and the initial condition G(x, y; t) — for t < 0. The Fourier transform 

~ roo 

G(x,y;cu)= dtG(x,y;t)e iuJt (2.3) 
Jo 

satisfies 

D(u)G = [-co 2 -d 2 x + V(x)] G(x, y; u) = 5(x - y) (2.4) 

and is analytic in the upper half uj plane. 

Waves propagating on curved space and observed by a distant observer should satisfy 
the outgoing wave condition. In terms of this, there are two classes of systems, 
(a) For non-singular systems, e.g., an oscillating relativistic star generating gravitational 
waves, the variable x, as defined by (1.5), is the radial variable of a three-dimensional 
problem, < x < oo. The left boundary condition is taken to be 

4>(x = 0, u) = , (2.5a) 



while the right boundary condition is one of outgoing waves 

<j>(x, u) oc e iuJX , a; ^+00 , (2.5b) 
where is the wavefunction in the frequency domain. 

(b) For the case of black holes, x given by (1.5) is a full line variable, — 00 < x < 00, where 
the left boundary x — > —00 is the event horizon. We have outgoing waves at both ends, i.e., 

<j>(x, u) oc e~ iujx , , (2.6) 

together with (2.5b). 

For simplicity, we shall write the formalism in terms of the half-line problem, and indicate 
briefly the simple changes necessary for the full-line case. We shall assume V — > as 
\x\ — > 00. 

Define two auxiliary functions f(u, x) and g(u, x) as solutions to the homogeneous equa- 
tion D(u)f(u,x) = D(u)g((j,x) = 0, where / satisfies the left boundary condition and g 
satisfies the right boundary condition. To be definite, we adopt the following normalization 
conventions: 

]im[e-^(^,x)] = l; (2.7) 

and for the half-line problem f(u,x = 0) = 0,f f (u,x = 0) = 1 while for the full-line 
problem \im x ^_ 00 [e' iUJX f (u , x)] = 1. Let the Wronskian be W{uS) = W(g, f) = gd x f — fd x g, 
which is independent of x. It is then straightforward to see that 

W{uj) ' X<y 
G(x,y;u) = < . (2.8) 

Wp) ' y<x 

The essence of the problem is, in one guise or another, how to integrate / from the left and g 
from the right, until they can be compared at a common value of x, at which the Wronskian 
can then be evaluated. The above choice of normalization is for convenience only, since a 
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change of normalization, / — > Af and g — > Bg, would result in W — > AB W with no effect 
on G. 

Consider the inverse transform of (2.3), and for t > attempt to close the contour in uj 
by a semicircle of radius C in the lower half-plane, where eventually one anticipates taking 
C — > oo. One is then led to consider the singularities of G for Im < (Fig. 2), from which 
different contributions to G can be identified as follows. 

(a) Tail contribution If the potential V(x) has a finite range, say V(x) = for x > a, 
then one can impose the condition (2.7) at x = a + , and integrate the equation through a 
finite distance to obtain g(u,x) for all x < a. Such an integration over a finite range using 
an ordinary differential equation cannot lead to any singularity in uj [16j . It is not surprising 



that if V(x) does not have finite support, but nevertheless vanishes sufficiently rapidly as 
x — > +oo, g(uj,x) is still regular. The necessary condition is that V(x) must vanish faster 
than any exponential |17|, and this is shown in Section IV. However, if V(x) has a tail that 



does not decrease rapidly enough, then g(u,x) will have singularities in uj. A tail that is 
exponential at large x, V(x) ~ e~ Xx , would in general lead to singularities in g(uj,x) for 
— Im uj > A/2. In so far as these singularities stay away from uj = 0, they have no effect on 
the late time behavior. However, if V(x — > oo) ~ x~ a (logx)P {f3 = 0, 1), then g(uj,x) will 
have singularities that lie on the — Im uj axis, and take the form of a branch cut going all 



the way to the origin ||17|| , as in the Schwarzschild case [|13|]. This contribution to G will be 
denoted as G^. It is the purpose of this paper to evaluate Gl for a broad class of V(x). 

For the half-line problem, f(uj, x) is integrated from x = through a finite distance, and 
hence does not have any singularities in uj. For the full-line problem, f(uj,x) is dealt with 
in the same manner. In all cases of physical interest, V(x) vanishes on the left either faster 
than an exponential, or precisely as an exponential [e.g. Schwarzschild spacetime under the 
transformation (1.5)]. For the former, f(uj,x) has no singularities in uj, while for the latter, 
there will be a series of poles, but at a finite distance from uj = 0. In both cases, there is no 
contribution to the late time behavior. 

(b) Quasinormal modes Secondly, the Wronskian W(uj) may have zeros at uj = ujj,j = 
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±1, ±2, . . ., on the complex to plane ||18|| . At these frequencies, / and g are linearly depen- 
dent, i.e., one can find a solution that satisfies both the left and the right boundary condition. 
Such solutions are by definition quasinormal modes (QNM's). The collective contribution 
of the QNM's is denoted as GQ(x,y;t). In the case of a Schwarzschild black hole, Gq 
has been evaluated semi-analytically by using the phase- integral method [T9|]. Since each 



Im ujj < 0, Gq decays exponentially as t — > oo, and is, therefore, irrelevant for the late time 
behavior. In some cases, GQ(x,y;t) is the sole contribution to the Green's function; this 
would of course be important, since the quasinormal modes would then be complete. This 
case will be discussed elsewhere HITf . 

(c) Prompt contribution Finally there is the contribution from the large semicircle at 
\u\ = C,C — > oo. Since \u\ is large, this term contributes to the short time or prompt 
response, and the corresponding part of the Green function will be denoted as Gp(x,y;t). 
This term can be shown to be zero beyond a certain time; thus it does not affect the late 
time behavior of the system. 

The three contributions Gl,Gq and Gp can be understood heuristically using a space- 
time diagram of wave propagation, as illustrated in Fig. 3, in which free propagation is 
represented by 45° lines, and the action of V(x) is represented by scattering vertices. The 
initial wavefunction is assumed to have compact support, so the propagation starts at some 
initial y, and one imagines observation at some x. Rays la and lb show "direct" propagation 
without scattering. These rays arrive promptly at x, and correspond to Gp. The ray labelled 
as 2 suffers repeated scatterings at finite x. For "correct" frequencies, the multiple reflections 
add coherently, so that retention of the wave is maximum. Moreover, the amplitude decreases 
in a geometric progression with the number of scatterings, which is in turn proportional to 
the time t. Thus (j) decreases exponentially in t; these contributions correspond to Gq. 

Lastly, for a wave from a source point y reaching a distant observer at x, there could be 
reflections at very large x', such as indicated by ray 3. The ray first propagates to a point 
x' ^> x, is scattered by V(x'), and returns to the observation point x, arriving at a time 
t ~ (x' — y) + (x' — x) ~ 2x'. If the potential V(x) ~ x~ a (logx) 13 with a > 2, then the 
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scattering amplitude goes as V(x') ~ V(t/2) ~ t _a (logt) /3 . Thus this contribution leads to 
Gl- In Section V, we shall see that this heuristic argument does not work for a = 2. In 
particular, a centrifugal barrier of angular momentum I corresponds to free propagation in 
three-dimensional space, and does not contribute to the late time tail. 

This heuristic picture requires two modifications, which we shall demonstrate both nu- 
merically and analytically in the following Sections. First, for V(x) = 1(1 + l)/x 2 + V(x) as 
x — > oo, where V(x) ~ x~ a (\ogx) 13 , the late time tail is due to V(x) and turns out to be 
*-(«+*) (log generically. The suppression by a factor t 21 , at least in the case a = 3, is 
known for specific black hole geometries @||. Second, if a is an odd integer less than 21 + 3, 
the leading term in the late time tail vanishes. For (3 = 0, the next leading term is expected 
to be £-( 2/ + 2a -2) ? while for (3 = 1, the next leading term is t _ ( 2/+Q ) without a logt factor. 

A minor technical complication should be mentioned at the outset. We close the contour 
by a semicircle of radius C (or other large contour). For any finite C, one can define the 
contributions from (a) the cut, (b) the poles and (c) the large circle. When C — > oo, the 
sum of these must tend to the finite limit G(x, y; t); however, there is no a priori guarantee 
that each of the three terms would individually converge to some limit. It can be shown 
that provided a regulation is introduced, then each term would individually approach a 
finite limit, and hence the tail, QNM and prompt contributions would each be well defined. 
The regulator could be e _arr , r —>■ + (for the part coming from Re u > 0). The need for a 
regulator comes from the region |Im u\ = 0(log |Re u\), |Re u\ — > oo, and has no bearing 
on the tip of the cut, near the origin, which controls the late time behavior. Therefore, the 
regulating factor will not be written out in the rest of this paper. 

To study the late time behavior, it then remains to study the spatial asymptotics of 
the potential and the consequent singularities of g. Before carrying out these analyses for 
various potentials in Sections IV to VI, we first explore the late time behavior numerically. 
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III. NUMERICAL CALCULATIONS 



In this Section, we present numerical investigations of the late time behavior. Eq. (1.1) is 
integrated numerically in time using standard techniques. Instead of the usual second-order 
differencing for the spatial derivative term, which could lead to erroneous results on account 
of a subtlety to be discussed in Appendix A, we have used a higher-order differencing scheme. 

We evolve 4>(x,t) for the half line x G [0,oo), with the boundary conditions = at 
x = and outgoing waves for x — > oo. The full line case [x G (— oo, oo), with outgoing wave 
boundary conditions for \x\ — > oo] is basically the same. From (2.1), the generic leading late 
time behavior is given by the integral / dyG(x, y; t)d t (f)(y, 0), so we take <ft(y, t = 0) = and 
use gaussian initial data: d t (f)(y, t = 0) = e~( y_2/ °) ^ . 

We study a broad class of potentials: 

V(x) = ^±V+V(x) , (3.1) 

Ob 

with 

,1/3 



V(x) = 



2T 







log 


(-)l 

\X J J 



x > x c , (3.2) 



for some a > 2, and some x and x c , which includes (i) power-law potentials {(5 = 0) and (ii) 
logarithmic potentials ((3 = 1). We study both integral and non-integral v. The logarithmic 
potential is interesting as it includes the Schwarzschild metric as a special case. For x < x c , 
we round off the potential by taking 



2u(u+l) 2x a - 2 f, fx + x r \^P 



x 2 + x 2 c ) (x a + a£) 



log 



2x 



o 



x < x c . (3.3) 



We have checked that the form of the late time behavior does not depend on the round-off. 
For all the numerical calculations shown below, we take x c = 100, y Q = 500 and rj = 50. 

We first consider power-law potentials V with v — and study the dependence of the 
late time tails on the parameters a and x Q . Fig. 4 shows log \<ft(xi,t) \ versus logt at a fixed 
point xi > x c for different values of a and x Q . The solid lines are numerical results while 
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the dashed lines are analytical results, which will be derived in the following Sections. For 
all cases studied, the tails are of the form t~ M with /i ~ a. For a given a, the magnitude of 
the tails is larger when x Q is larger. 

Next, we look at the case when V(x) is zero and investigate how the late time tail 
depends on v. It turns out that the result depends significantly on whether v is an integer. 
For non-integral z/, the tail is again of the form t~ M but with /z ~ 2 + 2v. On the other hand, 
when v is an integer, we see only quasinormal ringing and no power-law late time tail. This 
qualitative difference in the tail is observed even when v is altered by a small value, and 
this is demonstrated in Fig. 5, which shows log t)| versus logt for v — 1 and 1.01. 
This result suggests that the magnitude of the tail is probably proportional to sin z/7r, and 
therefore vanishes when v is an integer. 

Then, we study power-law and logarithmic potentials with nonzero integral values of u, 

i.e., v — I — 1,2, Two interesting results are found. First, logarithmic potentials often 

lead to logarithmic tails of the form t~ M logt, with /i ~ 21 + a. Thus, the late time tail 
is not necessarily an inverse power of t. To exhibit the logarithmic factor, we plot \<p\t 2l+a 
versus logt for several logarithmic potentials with / = 1 in Fig. 6. The existence of a logt 
factor is indicated clearly by the sloping straight lines at late times [see below for case (iii), 
which is an exception]. However, when we vary the parameter a continuously from 2.9 to 
3.1, the behavior of the late time tail changes discontinuously. This is seen from case (iii) in 
Fig. 6. Case (iii) is asymptotically flat, showing that there is no logt factor, and that the tail 
becomes a simple power law ~ £-(2Z+a)_ g uc h discontinuous jumps are observed whenever a 
goes through an odd integer less than 21 + 3. Interestingly, the Schwarzchild metric is such 
an exception, and the late time tail decays as t~( 2l+3 \ with no logt factor, as is well known. 

Discontinuous jumps are also observed for the power-law potentials at these values of 
a, and the change in behavior is more pronounced. The generic late time behavior for the 
power-law potentials is t~ M with /i ~ 21 + a, but when a goes through an odd integer less 
than 21 + 3, fi jumps to a larger value ~ 21 + 2a — 2. One such discontinuous jump is 
demonstrated in Fig. 7 with / = 1 and a = 2.9, 3.0 and 3.1. 
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These results and other cases studied but not shown here are summarized in Table 1. In 
the following sections, we present an analytic treatment for the various different potentials. 
The results so obtained will be shown to agree very well with the numerical evolutions 
plotted here. 

IV. POTENTIAL DECREASING FASTER THAN X~ 2 

The differential equation D(u)g(u,x) = 0, together with the boundary condition at 
x = +00, leads to the exact integral equation 

g(u,x) =e iwx - I™ dx' Smu{ - X ~ X '\ {x')g{LU,x') . (4.1) 

Jx LO 

In Born approximation, the factor g(u, x') in the integrand may be replaced by e lux> , so (4.1) 
can be reduced to 

g{u,x)=e iu *-I{u,x) , (4.2) 

where 

I(u,x)= / dx -V(x')e lu}X . (4.3) 

Jx LO 

Note that if V(x) is of finite range, then there can be no singularities in u. 

Suppose the potential decreases like an exponential, i.e., V(x) oc e™ Ax , then we have 

e (—\+iu))x 

implying a pole at u = —i A/2 on the - Imw axis. Heuristically, if the potential decreases 
faster than any exponential, then A — > 00 and the pole effectively disappears, resulting in 
an analytic g(u,x). As a result, there is no tail contribution to the Green's function. For 
the exponential potential, g(u, x) can actually be solved exactly without using the Born 
approximation ]I7[]. The pole structure and the outgoing wave condition can then be seen 
explicitly |20(| . In any event, since these singularities stay away from u = 0, they will not 
contribute to the late time behavior. 
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Next, we consider potentials decreasing less rapidly than an exponential, and in particular 
v — and a power-law V, viz., 



n.<-) = -A- 



X* \ X 



(4.5) 



We take a > 2 and the coefficient K will be absorbed into the scale parameter x Q . When a 
equals an integer n, we have (Appendix B) 

,({2iux )^ 



I(uj, x) 



■ ) \ 

I (n-l)! 
1 



n-2 y 

7 — E h log(— 2iuix) 



+ 



x, 



m=l 
n-2 oo 



E 



(2ia;a;) ? 



(4.6) 



n-\\x) m=0 ^ n _ 2 m!(m + 2-n) j 

where 7 is Euler's constant. The appearance of e~ luJX in J and thus in g(u,x) indicates 
that a right-propagating wave at x — +00, when extrapolated to finite x, contains a small 
admixture of left-propagating wave. With reference to ray 3 in Fig. 3, it is clear that this is 
the term responsible for the late time behavior. The logarithm in (4.6) causes a branch cut 
along the - Imw axis (see Fig. 2). To describe the cut, let 



g±{-i(T,x) 



lim g(—ia ± e, x) 



for a real and positive, and 



A(a,x)= g+(-ia,x) -g_(-ia,x) 

= — \im + [I(—ia + e, x) — I(—ia — e, x)] 

Then from (4.2) and (4.6), the strength of the cut is 



(4.7) 



A(<7, x) ~ 2ni 



.(2<7X ) 



n-2 



(n-1)! 

For non-integral a, we have (Appendix B) 



(4.8) 



I(cu, x) 



(a-l) 



-2iujx ) a - 2 T(2 -a)- 



x 



0.-2 00 



E 



{2iujx) 7 



(4.9) 



xj ^ Q m\(m + 2- a) 

In this case, the term (— 2iu)x ) a ~ 2 causes a cut, again along the - Imw axis. The strength 
of the cut is 
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A(a,x)^2ni {2 ™° )a ~ 2 e-™ , (4.10) 
r(a) 

which reduces to (4.8) when a — > n. In this sense, the result for non-integral a is more 
general. 

The Green's function G(x, y; uo) as given by (2.8) will likewise have a cut along the -Imw 
axis. The values on either side are given by 

G *(«-*«> = w i9± ,f) • (411) 

where henceforth we take y < x. A little manipulation gives the discontinuity 

AG(x, y; a) = G+(x, y; u) - G-(x, y; u) | w= _ i(7 
W{g-,g+) 



W(g + ,f)W(g_,f) 

From (4.2) and (4.9) 



(4.12) 



W^U-, ^+) U=-«t r( x • (4-13) 

Since W^((?_,(? + ) is independent of x, in order to derive (4.13) it is only necessary to know 
g± at one value of x. For a potential decreasing faster than x~ 2 , the Born approximation 
is valid for very large x, and we have used the resultant (4.9) only in that region. Nowhere 
else do we rely on the Born approximation. Therefore, these results are exact, and do not 
suffer from any inaccuracies arising from the use of the first Born approximation. In other 
words, the t — > oo behavior relates only to scattering at x' — > oo, for which the first Born 
approximation is accurate, provided a > 2. It remains to evaluate W(g±, f). To be specific, 
consider the half-line problem and evaluate the Wronskian at x — 0, 

W(g±J)=g±(u,x = 0). (4.14) 
Substituting (4.13) and (4.14) in (4.12), we have 

AG(x, y; a) * J ^^Tl f(-i*,x) f(-i*,y) 
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The contribution arising from the cut along the — Im 00 axis is 

G L (x, y; t) = -i —AG(x, y; -ia^ . (4.16) 
Jo 2n 

The t — > 00 behavior is controlled by the a — > behavior of (4.15). A little arithmetic then 
leads to 

GAx y t) ~ (4 17) 

Let the initial conditions be <j)(y,0) = <f>o(y), d t (f)(y,0) = (f>i(y), and for simplicity we 
assume O , 4>i to have compact support, or at least to vanish at infinity faster than any 
exponential. Then the late time behavior is 



na— l„a— 2 r 

#M)- ^ 



/ dyf(0,y)Mv) 
Jo 



Sf ■ (418) 



The term due to <p is proportional to d t G; hence it goes as t~( a+1 \ and becomes subdominant 
when t — > oo. 

To compare the theoretical result (4.18) with the numerical evolutions in Section III, 
we obtain the constant <?(0, 0) and the spatial function f(0,x) by integrating the time- 
independent equation D(u — > 0)/ = D(lu — > 0)g = 0. The two results agree very well when 
logt is large (see Fig. 4). In Fig. 8, we plot versus a; at a time when the t~ a behavior 
is observed, showing excellent agreement also in the x-dependence of at late times. 



V. INVERSE SQUARE LAW POTENTIAL 

In this Section we consider inverse-square-law potentials (V = 0): 

V(x)= y 2 1 , x>x c . (5.1) 

Ob 

We assume the potential is repulsive and take v > 0, and define p = u+1/2. For V oc x~ 2 , the 
Born approximation is not accurate at any x (unless the coefficient of the potential is small), 
so the derivation in Section IV is not valid. However, in this case, an explicit solution can be 
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written down. The solution satisfying the outgoing wave boundary condition as x — > +00 
is exactly 



g{u,x)=e^nf-^Hy{ux). (5.2) 
Since H^{z) = J p {z) + iN p {z) and 

W- f ^T)(^ + ... , (5.3a) 

N p (z) ~ -^.fy-P + . . . , (5.3b) 

the origin is a branch point for non-integral v. The symmetry in the present problem dictates 
that the cut should be placed along the — Im ui axis, and some algebra leads to 



g + (-ia,x) = a /^e^ /2 [2 cos pnH^iiax) + e' ifm Hf\iax)] , (5.4a) 



g-(-ia,z) = ^^^l\e-^Hf{iax)\ , (5.4b) 

where the functions H^ ,2 \iax) have no discontinuity for a > 0. Then 

W(g_,g + ) | w =_ iff = -4msinz/7r , (5.5) 

which vanishes for integral z/. This is readily understood since V(x) = l(l + l)/x 2 represents 
a pure centrifugal barrier, which corresponds to free propagation in three-dimensional space 
and does not cause any scattering. Note that for v — > 0, the potential is weak and the 
Born approximation derived in Section IV should be valid. It can be readily checked that 
the v — > limit of (5.5) indeed agrees with the a — > 2 limit of (4.13). [For a — 2, it is no 
longer possible to absorb the magnitude K of the potential into the definition of x Q) so (4.13) 
should be multiplied by K = v{y + 1) ~ v in this case.] This agreement may be regarded 
as a check on the formalism in Section IV. Thus the dependence W(g-,g + ) oc cr n ^ 1 remains 
valid for n = 2, but the amplitude vanishes. 
For 00 = —ia — > 0, 
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g_(-ia,x)~ —e*-»l 2+1 I^T{p) 

V 71 " \ ^ 



2 



for x > x c . Since g± oc a; 



f'(x c ) + -f(x e ) 



9±(x c ) , 



(5.6a) 
(5.6b) 

(5.7) 



and we have evaluated the Wronskian at x — x c . As discussed in Section II, f(u,x) is 
analytic in u, which implies f(u>,x) = Y^=o fn{x)oj n . From the Wronskian, it can be seen 
that f (x) ^ 0. Thus, 



hm W(g + ,f)W(g-,f) = - 



f'o(x c ) + —f (x c ) 

Xr 



2 T( P ) 2 (<jg\- 2v 

71 V 2 ) 



Using (4.12), we have 
AG(x,y;a) 



Am sin vix 



[fo( X c) + (v/x c )Mx c ) 

and the late time behavior is 

2 sin vix 



x r \ 2v a 2u+1 



2 ^ 2 J r(z/ + |) 2 



T^fo(x)fo(y) 



G L (x,y;t) ~ 



r(2 ^ + 2) ^fr<^>/ (*)/„(y) 



(5.8) 



(5.9) 



(5.10) 



Two features are worthy of remark. First, whereas the lowest order Born approximation 
in Section IV would have predicted a time-dependence of t~ 2 for a potential V ~ x~ 2 , 
the power is in fact £-( 2 + 2l/ > for non-integral v. For v <C 1, the two agree, whereas for 
z/ = 0(1), the difference arises from higher-order Born approximation. Secondly, the power- 
law dependence in t vanishes for integral v and the late time behavior is then dominated by 
the exponentially decaying QNM's. 

The late time behavior is 

2sin(z/7r) T(2v + 2) 



(f)(x,t) 



fo(x c ) + vfo(x c )/x t 



r|(y) 2 ^" (2+2l/) /o(x)/ o °°^/o(y)0i(y) . (5.11) 



This analytical result for <p{x,t) is evaluated and plotted as the dashed line for v = 1.01 in 
Fig. 5. The dashed line is indistinguishable from the solid line (obtained from numerical 
evolution) when logt > 8. On the other hand, for v — 1, <f> decreases much faster than 
t~( 2+2u \ in fact faster than any power law, as derived above. 



19 



VI. COMPOSITE POTENTIAL 



We now consider potentials of the following form: 

V ( x ) = l l±±l + V(x), x>x c , 1 = 1,2,... (6.1) 

x z 

with V(x) given in (3.2). We first take a > 2 to be non-integer, and results for integral a 
will be obtained by taking the limit a — > n, where n is a integer. For I = and j3 = 0, V(x) 
reduces to the case discussed in Section IV. 

We see from Section V that if V is treated as a perturbation, the lowest order Born 
approximation no longer gives correct results. Rather, we should consider V(x) as the 
perturbation. The solution for the inverse square potential, 1(1 + l)/x 2 , is given exactly by 

g {0) (uj,x) =i l+ \uox)hf\ujx) , (6.2) 

where is the spherical Hankel function of the first kind. Applying the Born approxima- 
tion now gives, in analogy to (4.2), 

9{u,x) ^ g (0) (uj,x) -I(u,x) , (6.3) 

where 

I{u,x)= dx'M(x,x';u)V(x')g {0 \u,x') , (6.4) 

Jx 

with the zeroth-order Green's function being 

M(x,x';u) = l -^-xx' [hf\ux')hf\ux) - hf\ux)hf\ux')] , (6.5) 

and fif* is the spherical Hankel function of the second kind. 
First consider a power-law V. After some algebra, we get 



i l ux 



I\{uj, x)h\ 1 \ujx) — hiu, x)h\ 2 \ujx) , (6.6) 



where 
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h{u,x) = (ux ) a - 2 Hdu hl {u l h [ {U) , (6.7) 



/ uix U 

r°° h^(u) 2 

,n-2 / ./.. "/ K U ) 
II' • 2 



I 2 (u,x) = (Lux ) a - 2 / du-^-f- . (6.8) 
While Ii is straightforwardly evaluated to be 

h(u,x) = 2- 2l (ux)- 21 - 1 (?°) Yl a m (l,a)(cux) 2m , (6.9) 



x a-2 i 



^ X ' m =0 

I2 is somewhat more difficult to obtain. We find 



I 2 (u,x) = (-l) l+l 2iC{l, a) 7 2 ^ {-2iux ) a - 2 

(a - 1) 

\ a-2 00 



+(-i)Wr 2 ^ - EMUIH" , (6.10) 



m=0 



where 

^.-)=n^^. '=1,2,... , (6.11) 

and C(Z,a) = 1 for / = 0. [See Appendix C for the derivation of (6.9) and (6.10) and the 
definitions of a m and b m ] The term (— 2iux ) a ~ 2 causes a cut in g(u,x), which lies along 
the negative imaginary axis, with a strength 

A(a,x) -n(-ty + \2ax ) a - 2 ^^-(2ax)h < i 2 \-tax) . (6.12) 

r(a) 

Thus the factor W(g-, g + ) is given by 

W(g_,g + )l = _ ia ~ W(g^\ A) = -4f W (2 ra /- 2 ?^ . (6.13) 

r(«) 

In the case of a being an integer n > 3, we get 

W%-,<? + ) = -47ra(2^- 2 ^^ . (6.14) 
Approximate W(g±, f) by W(g(°\ f), and we have 



\\mW(g + J)W(g^f) = a- 2l g 2 , 

a — >U 
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where g Q = \\m^^ Q [{iu}) l W {g° , /)] is a non-vanishing constant and reduces to <?(0, 0) for 
Z = 0. Thus, 



AG(,, V ; -ic) g - Ti 2»,;-V a+ »-' ^ f /H<7 ' ^{'"'^ y) , (6.15) 



and as a result, 



9o r («) 
f(0,x)f(0,y) F(a) 

= ~2 CU")-^ ( 6 - 16 ) 

for both integral and non- integral a. 

We note the interesting result that C(l,a) = when a equals an odd integer less than 
21 + 3, in which case the late time tail vanishes in first Born approximation, and higher order 
corrections have to be considered. This implies multiple scatterings from asymptotically far 
regions are important in such cases. We shall not go into the details of the higher order 
calculations, but from dimensional analysis, we expect 



G L (x,y;t)~ t -^ 



+*(f) + 



(6.17) 



where A{ 's are constants proportional to the i -order scattering amplitude. 

Generically the late time behavior arising from the first Born approximation is linear in 
the potential. The potential with a logarithmic V is obtained by taking (— d/da) on that 
with a power-law V . So we immediately obtain from (6.16) that for logarithmic V 

Gl(x, y; t) = fa, a)F(a) 1 ^ - A [C(l, a)F(a)] ^} /( °' ^ y) . (6.18) 

Hence, the leading terms are t~^ 2l+a \c\ogt + d), except that c oc C(l, a) vanishes when a is 
an odd integer less than 21 + 3. 

Analytical results obtained from (6.16) and (6.18) with no adjustable parameters [except 
in case (ii) in Fig. 7; see below] are plotted as dashed lines in Figs. 6 and 7. The agreement 
between numerical evolutions and analytical results is perfect. For case (iii) in Fig. 6, the 
vanishing of the leading term implies that the asymptotic slope should be zero (i.e., no 
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log t, but only a pure power, whose magnitude is determined), and this indeed agrees with 
the numerical results, with no adjustable parameters. For case (ii) in Fig. 7, the leading 
term vanishes, and the dashed line represents the next leading term arising from multiple 
scatterings, whose time dependence is determined, but whose magnitude has been left as an 
adjustable normalization. In Fig. 9, we show the dependence of on x at late times; perfect 
agreement is again found. 

VII. CONCLUSION 

In this paper we have studied the dynamical evolution of waves on a curved background 
that can be described by the Klein-Gordon equation (1.1) with an effective potential. This 
includes linearized scalar, electromagnetic and gravitational waves on a Schwarzschild back- 
ground as special cases. Our work extends, places in context and provides understanding 
for the known results for the Schwarzschild spacetime. 

We have given a systematic treatment of waves described by (1.1) using a Green's func- 
tion formulation. The Green's function consists of three distinct components, Gl, Gq, 
and Gp, coming from the three contributions to the contour integral in frequency space. 
Each component leads to a distinct wave phenomenon in physical spacetime. Gq leads to 
quasinormal ringing, Gp leads to prompt responses, and Gl leads to the late time tail phe- 
nomenon. Fig. 1 shows clearly these three contributions, while Fig. 3 provides a heuristic 
understanding in terms of wave propagation in a spacetime diagram. 

We further focus on the late time tail phenomenon. Both numerical and analytic results, 
which agree perfectly with each other, are presented. For a broad class of potentials, we have 
demonstrated that the tail is related to the Green's function having a branch cut along the 
- Imw axis, with the asymptotic late time behavior controlled by the tip of the cut. Using 
a Born analysis, we determine the strength of the cut in terms of the spatial asymptotics 
of the potential. This gives analytic understanding of both the magnitude and the time 
dependence of the late time tail, not just for the Schwarzschild case, but also for a wider 



23 



class of models. We have applied our formulation to a few interesting cases, including the 
inverse- square-law potential, other power-law potentials decaying faster than inverse square, 
potentials containing a logarithmic term, and composite potentials with a centrifugal barrier. 

We find that the well-known late time power law behavior of waves on a Schwarzschild 
spacetime turns out to be an exceptional case. In general, for potentials going as a centrifugal 
barrier of angular momentum / plus V(x) ~ x~ a logx, the generic late time behavior is 
~ £-( 2 '+ a ) logt, which is not a simple power law. However, when a is an odd integer less 
than 21 + 3, as in the case of Schwarzschild, the coefficient of the leading term vanishes, and 
the next leading term is £-( 2/ + a ) without a logt factor. This discontinuous dependence of 
the time dependence on a is verified numerically and shown in Fig. 6. 

There are various implications, colloraries and interesting future extensions discussed 
throughout the paper. Appendix A outlines a subtlety in the numerical calculations of the 
late time tail behavior (or more generally speaking the accurate time dependence of solutions 
of partial differential equations) in the timelike direction. 
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APPENDIX A: GHOST POTENTIAL 



In this Appendix we consider a subtlety in the numerical solution of (1.1), which we 
write as 

d 2 <f> — L [(/>] = d 2 x J> - Vcf> . (Al) 

The numerical solution does not satisfy (Al), but instead satisfies a modified equation with 
L replaced by 

L\^=D 2 X <\>-V<\>, (A2) 

where D 2 X is a finite-difference approximation to d 2 , e.g., in a second-order scheme, 

[<p(x + Ax, t) - 2(f)(x, t) + (f)(x - Ax, t)} 



1 

12' 



- d 2 x (j)+—(Ax) 2 dt(j)+--- . (A3) 



For the class of problems at hand, the solution takes the form 

<f>(x,t) = r»f(0,x) + --- , (A4) 

where \i is some exponent, f(cu,x) is the solution regular at the origin, and the omitted terms 
are subasymptotic as t — > oo at fixed x; see e.g. (6.16). (The convergence is not uniform in 
x, and moreover, a logt factor can be incorporated without changing the argument below.) 
To lowest order in Ax, this satisfies (Al), and noticing that d 2 <j) is down by a factor t~ 2 , we 
see that d 2 (p ~ V<j) + • • -. This allows us to eliminate all spatial derivatives on higher than 
second order, and a little arithmetic leads to 

L[4>] — L [0] + l(Ax) 2 [(V" + V 2 )<P + 2Vy] + • • • , (A5) 

where ' denotes d x . 

Two further simplifications are possible. First, because of (A4), the d 2 (j) term is sub- 
asymptotic as t — > oo, and consequently the error caused by finite differencing in t can be 
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ignored, at least to order (At) 2 . Secondly, assume that the potential is dominated at large 
x by the centrifugal barrier, then 

f(0,x)~Kx l+1 (A6) 

[The function /(0, x) corresponds to the radial solution to Laplace's equation, and is regular 
at the origin. In general, it contains both the decreasing and the increasing component. At 
large x, the latter dominates, and (A6) then follows. The behavior (A6) has been verified 
numerically] As a result, <p' ~ [(/ + l)/x](f>, and it is possible to write the extra terms in 
(A5) purely in terms of without 0': 

Z[0] ~ L[<i>]+V gh (x)<j> + -- - , (A7) 

where the ghost potential is 

oc (Ax) 2 x" 4 . (A8) 



V gh (x) = ^(Ax) 2 



v" + v 2 + 2 ( Z + 1 V 

x 



Thus if V vanishes at infinity faster than x~ 4 , say as x~ a (a > 4), then instead of seeing 
the true late time behavior t~( 2l+a \ the numerical solution based on second-order spatial 
differencing would display a late time behavior £~( 2?+4 ) due to the ghost potential. It is easy 
to be misled by the ghost behavior unless one has analytic expectations to compare with. 

As an even clearer example, consider a pure centrifugal barrier, say with 1 = 1. The late 
time behavior should be exponential (since there should be no cuts along the - Imw axis, 
and the leading large t contribution comes from the quasinormal modes). The top line in 
Fig. 10 shows the result of a numerical calculation with second-order spatial differencing; the 
straight line portion gives a dependence t~ 6 , which is the ghost behavior as explained above. 
The lower line in the same figure shows the result when Ax is halved; the t dependence 
remains unchanged, but the magnitude is decreased, roughly by a factor of 4, as expected. 
The change with Ax shows that this tail is a numerical artifact. 

In short, k-th order spatial differencing leads to 

V gh {x) oc (Ax) fc x~( 2+fe ) , (A9) 
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which could mask the effect of V if the latter decreases as x — > oo more rapidly. This problem 
could have led to errors, or at least ambiguities, in some numerical experiments. In all the 
numerical results presented in this paper, we have taken care to use spatial differencing of a 
sufficiently high order that V g h will not dominate over V, and have moreover checked that 
the t — > oo asymptotics claimed are stable against (a) decreasing Ax, and (b) increasing the 
order of spatial differencing. 



APPENDIX B: 



For a power-law potential of the form (4.5) with K — 1, using (4.3), we have 



I{uj, x) 

which is easily shown to be: 



00 f sinu(x — x') x a 2 



dx' 



o £—lU)X 

x a 



I(u>, x) 





r°° du 


^iujx 


r°° due 2lu 




2i \ 


Jwx U a 




Juix U a 





After evaluating the first integral and integrating the second one by parts, we get 

J-2icux o r^ 



I(uj, x) = —e 



(a-1) 



-J(U!X) 



where 



J(ujx) = 



OO g 

du 

-2iujx W 



a-1 



If a is an integer n, we have [21]: 



J{ujx) 



:-i) 



n-l 



(n-2)! 



ra-2 



7 — h log(— 2iuix) 



m=l 



m 



i \n— 1 oo 
m=0 



(2iujx) 7 



[B2) 



(B3) 



(2iux) n 2 ^_ n f^„_r, ml(m + 2 — n) 



(B4) 



where 7 is the Euler's constant. 



On the other hand, if a is non-integer, we have J21 



J(ujx) = T(2 - a) + 



E 



{2iujx) 



(2i(jjx) a - 2 ^ m\(m + 2-a) 



(B5) 



Substituting (B4) and (B5) into (B3) then gives us (4.6) and (4.9) respectively. 
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APPENDIX C: 



The spherical Hankel functions of the first and second kinds can be written as f22|| : 

hf\u) = r*-w» £ {l u t m \ (-2iu)- m , (ci) 

m =o m\{l-m)\ 



and 

hf\u) = i l+x u x e- iu V ^± m)! . , (2zn)- m . (C2) 
Using (CI) and (C2), Ii as defined in (6.9) is easily evaluated to be: 



h(uj,x) 

2~ 2; (Xo\ a - 2 X [(-1)"* + (-1)""] (2/ - mQ!(2/ - m 2 )! +TO2 

(ws) 2W U/ mi ^ =0 2(2/ + «-l -m x -m 2 )(/-m 1 )!(/-m 2 )!m 1 !m 2 ! 1 J 



Now since (— l) mi + (— l) m ' 2 vanishes when mi and m 2 are not both odd or both even, we 
have only even powers of ux in the double sum. Hence, we can define the double sum as 

i i (2i) mi+m2 [(-l) mi + (-l) m2 ]2(2/-m 1 )!(2/-m 2 )! 



]T J2 — — — — — 



mi=0 m,2=0 
I 

\2m 



mi+ni2 

(21 + a — 1 — mi — m 2 )(Z — mi)!(I — m 2 )!mi!m 2 ! 



LUX) 



■ 



m=0 

and get (6.9). 

To evaluate 7 2 , we consider more generally the following integral: 

J(/3, z) = ^ ^-dt , (C3) 

where v(t) is an analytic function and lim^oo v(t)t n = for any positive number n. Assume 
that v(t = 0) 0, so that the integrand is singular at t = for f3 > 0. Let the Taylor 
expansion of v(t) be J2^=o c nt n , which converges for t G [0, oo). 
Rewrite X((3, z) for (3 > 1 as follows: 



z(M=/ vy & dt + L -=y- , (C4) 
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where N is the largest positive integer less than (3. It is then straightforward to show that 

/•oo v (f\ _ ,n oo „n+l-/3 

m z) = i m ^L it _ £ ^_ , (05) 

Note that the integrand on RHS is well behaved at t — 0. We shall show that this integral 
is, in fact, equal to 

v(t) 



where the method of analytic continuation has applied to define it for j3 > 1. 

First, by considering the muliti- valued property of the function t a across the branch cut 
going from t = to t = oo, one can argue that for noninteger f3 

MP) = Tp^—TT I V 3dt ■ (C6) 
2i sin pir Jc W 

The integral goes along a contour C enclosing the positive real t-axis in the clockwise direc- 
tion as shown in Fig. 11. This representation has the advantage of being analytic in (3 for 
both positive and negative values of (5 and is readily applicable as an analytic continuation 
of the original integral. 

Second, by using a similar trick, one can also prove that 

f°° vW-T^Cnt" e"f» f t^W-ESf/,, (r7) 

Jo ¥ dt -2u^Jc tP ^ ■ (C7) 

Moreover, the integral 



^ n= \ " dt 

C 



vanishes, which can be seen by deforming the contour C into a infinitely large circle enclosing 
the origin, as shown in Fig. 11. This then completes our proof for the equality 

v(t) - E^o 1 c n t n 



and hence 



From (6.8), 



= jf V{ ) \r° Cn dt , (C8) 



oo n +l_/3 

mz)=up)-E^Y—5 ■ (C9) 

n=0 M 
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h(co,x) = (ujx ) a 2 J 2 (lux) , 



(CIO) 



where 



J 2 (cox) = / dt 



h?\tf 



Hx-2 



(Gil) 



We let 



w(t)= h?\t) 



.21+2 



(C12) 



so that w(t) is analytic in t. Using (CI), we have 

on U\— C — V* VV 1 \fci+fc 2 /o„-\n~fci-fc 2 (j + feQKj + ^2)!(~l) fcl+fc % n+2i-fc 1 -fc 2 

oo 

= (-l)w£d m t m . (C13) 

m=0 

Then 

J 2 (cx)= £ d t^ , (C14) 

which is of the form of (C3) with v(t) = w(t), (5 = 21 + a and z = uox. Hence using (C9), 
we get 

oo J ( \m+l-2l-a 

Mux) =i A + (-iy £ dm \7 - — , (cis) 

^ m + 1 - 21 - a 



Since X A is the analytic continuation for the case a < —21, it can be looked up in tables p2 

'+i(_^-i r (^)r(^) 2 ] 

~2"7r r(2 - a) 



Ja (-i^Hr 1 r(^)r(g^)^r(M) 



After some algebra, we get 

1 A = (-iy(-2ir^C(l, a)ffi^ (C17) 

(a-1) 

with C(l,a) defined in (6.11). Substitute (C15) and (C17) into (CIO) and after some sim- 
plification, we obtain (6.10) with 

b m = - : dm 9] • (C18) 

m + 1 — 21 — a 
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Table 1 



V(x), x — > oo 




4>{t), t — > oo 


(K/x 2 )(x /x) a 


all real a > 2 




v(v+ l)/x 2 


integer v 


e -*/T 


all non-integer v 


£-(2i/+2) 


l(l + l)/x 2 + (K/x 2 )(x /x) a 
I: integer 


odd integer a < 21 + 3 


f - * 4 , /i > 2/ + a 


all other real a 


£-(2i+a) 


1(1 + l)/x 2 + (K/x 2 )(x /x) a \og(x/x ) 
I: integer 


odd integer < 21 + 3 


£-(2J+a) 


all other real a 


r (2J+a) logt 



Table 1. Behavior of late time tails for various potentials 
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FIGURE CAPTIONS 



Fig. 1 Generic time dependence of (f>(xi,t) at a fixed spatial point x±, evolving from an 
initial gaussian and initial = 0. The time evolution is governed by the Klein- 
Gordon equation (1.1). The potential in this case is V(x) ~ 1(1 + l)/x 2 + (\ogx)/x 3 
with 1 = 1. (a) Shorter time scale to show the prompt contributions; (b) longer time 
scale to show the late time behavior. 

Fig. 2 Singularity structure of G(x, y; uj) in the lower-half cu-plane, and the various contri- 
butions to the Green's function G(x,y;t). 

Fig. 3 A spacetime diagram illustrating heuristically the different contributions Gl, Gq 
and Gp (see text for definition) to the Green's function. 

Fig. 4 log \(f>(xi, t)\ versus logt for potentials with v — and a power-law V, which go as 
x^ 2 /x a as x — > oo. (a) (i) x — 1, a — 3; (ii) x Q = 100, a = 3 and (b) (i) x Q = 1, 
a = 2.5; (ii) x Q — 1, a — 5.7. Solid lines are numerical results, which agree very well 
with the analytical results (dashed lines). 

Fig. 5 log \<f)(xi,t)\ versus logt for an inverse-square-law potential V(x) ~ v[y + l)/x 2 . (i) 
v = 1.01 and (ii) v — 1. Solid lines are numerical results. For z/ = 1.01, decreases as 
with /i ~ 4.03 which is in excellent agreement with the analytical result of t~( 2 + 2i/ ) 
(dashed line), derived in Section V. For v — 1, decreases faster than any power 
law, again agreeing with the theoretical prediction. For clarity, case (ii) is shifted 
downwards by 8.0. 

Fig. 6 A\<f)(xi,t)\t 2l+a versus logt for several potentials V(x) ~ 1(1 + l)/x 2 + \ogx/x a . (i) 
I — 0, a — 3; (ii) I = 1, a = 2.9; (iii) Z = 1, a — 3; (iv) I = 1, a = 3.1. For clarity, the 
data are multiplied by a constant A with (i) A = 10~ 9 ; (ii) and (iii) A = 5.6 x 10~ 10 ; 
(iv) A = 5.86 x 10~ 10 . The sloping straight lines observed at late times in cases (i), (ii) 
and (iii) indicate the existence of a log t factor. Case (iv) is an exception, in which the 
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nearly zero asymptotic slope shows that the tail is a simple power law. The numerical 
evolutions (solid lines) are indistinguishable from the analytical results (dashed lines) 
for \ogt > 9.5. 

Fig. 7 log \<f)(xi,t)\ versus logt for various potentials of the form 2/x 2 + l/x a as x — > oo. 

(i) a = 2.9; (ii) a = 3.0; (iii) a = 3.1. To make the three sets of lines stagger, vertical 
shifts have been applied. The generic late time behavior is t" M and /i ~ 2 + a except in 
case (ii) in which fi jumps discontinuously to 2a. Solid lines are numerical evolutions 
while dashed lines are analytical results, derived in Section VI. In cases (i) and (iii), 
the analytical results are completely specified both in form and magnitude; in case 

(ii) , only the form is known and the magnitude is fitted to the numerical result. 

Fig. 8 A(j)(x, t£) versus i at a fixed time £l for the same potentials plotted in Fig. 4. 
The numerical evolutions (solid lines) agree very well with the theoretical estimates 
(pluses), (i) x — 1, a — 3; (ii) x Q = 100, a = 3; (iii) x Q — 1, a — 2.5; (iv) x Q = 1, 
a = 5.7. The data are multiplied by a constant A for clarity, (i) A = 5 x 10 4 ; (ii) 
A = 10 6 ; (iii) A = 10 3 ; (iv) A = 2 x 10 15 . 

Fig. 9 A(p(x,t L ) versus x for various potentials of the form 2/x 2 + (cilogx + c 2 )/a; a as 
x — > oo. Solid lines are numerical results while the pluses are theoretical estimates. 

(i) Ci = 1, c 2 = 0, a — 3; (ii) c± — 0, c 2 = 1, a — 3.3; (iii) c± — 0, c 2 = 1, a — 4; (iv) 
Ci = 0, c 2 = 1, a — 5. The data are multiplied by a constant A for clarity, (i) A = 10 8 ; 

(ii) A = 10 10 ; (iii) A = 10 12 ; (iv) A = 10 15 . 

Fig. 10 Ghost behavior observed in the numerical evolution of (f>(xi,t) using a second- 
order spatial differencing scheme. The potential used in this case is V(x) ~ 2/x 2 as 
x — > oo and the expected late time behavior is an exponential decay as shown in Fig. 5 
(which uses a higher-order spatial differencing scheme) and Section V. In the present 
numerical evolution, however, a straight line portion indicating a time dependence of 
t~ 6 is observed. The magnitude of this "ghost" tail decreases by a factor of 4 when the 
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spatial interval Ax is halved (see Appendix A), showing that this tail is a numerical 
artifact. 

Fig. 11 Contour C for evaluating the integral Ta defined in (C6). 
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